Fluctuation theorem and large deviation function 
for a solvable model of a molecular motor 
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We study a discrete stochastic model of a molecular motor. This discrete model can be viewed 
as a minimal ratchet model. We extend our previous work on this model, by further investigating 
the constraints imposed by the Fluctuation Theorem on the operation of a molecular motor far 
from equilibrium. In this work, we show the connections between different formulations of the 
Fluctuation Theorem. One formulation concerns the generating function of the currents while 
another one concerns the corresponding large deviation function, which we have calculated exactly 
for this model. A third formulation concerns the ratio of the probability of observing a velocity i; to 
the same probability of observing a velocity —v. Finally, we show that all the formulations of the 
Fluctuation Theorem can be understood from the notion of entropy production. 

PACS numbers: 87.15.-v, 87.16. Nn, 05.40.-a, 05.70.Ln 



I. INTRODUCTION 

The living cell has evolved a diverse array of pro- 
teins, which can perform a variety of chemical functions. 
These proteins catalyze chemical reactions, control key 
processes like recognition or signalling, or act as molec- 
ular machines. Biological systems containing these ele- 
ments are typically described as active because they are 
in a non-equilibrium state as opposed to non-active sys- 
tems which can be considered at thermal equilibrium. 
Recently there has been a lot of interest in the non- 
equilibrium properties of active biological systems, such 
as the hair bundle of the ear [l|J3 j active membranes H , 
active gels [1] , active networks U, active lipid clusters [a| 
and living cells Q--- 

Clearly, one step towards understanding active biologi- 
cal systems, starts with the understanding of single active 
proteins, such as for instance molecular motor proteins. 
These molecular motor are nano-machines that convert 
chemical energy derived from the hydrolysis of ATP into 
mechanical work and motion [8, 9]. Important exam- 
ples include kinesin, myosin, RNA polymerase, and the 
Fq/Fi rotating motor. Mechano-transduction (i.e., the 
process of conversion of chemical energy into mechani- 
cal work) in motors has been described theoretically by 
ratchet-models (lol . [Tl| , which rely on the fruitful concept 
of broken symmetry. According to the Curie principle, 
directed motion requires to break the spatial symmetry 
and the time reversal symmetry associated with equilib- 
rium (detailed balance). The spatial symmetry is broken 
by the asymmetric interaction between the motor and 
the filament, while the time reversal symmetry is broken 
by chemical transitions, which break locally the detailed 
balance condition. From the continuous ratchet models 
described by Langevin equations ,10| . it is possible to 
construct discrete stochastic models of molecular motors 
by considering only localized discrete transitions as ex- 



plained in Ref. J^] . These discrete stochastic models are 
interesting because they are minimal, in the sense that 
they contain the main physical picture of ratchet mod- 
els while being more amenable to precise mathematical 
analysis [H, 111 [H, [3 . 

Recent advances in experimental techniques to probe 
the fluctuations of single motors pro vide wa ys to gain 
insight into their kinetic pathways [l7l [Tsl . [T9ll20|| . How- 
ever, a general description for fluctuations of systems 
driven out of equilibrium, and in particular of motors, 
is still lacking. Recently, the Fluctuation Theorem (FT) 
has emerged as a promising framework to characterize 
fluctuations in far-from-equilibrium regimes. This theo- 
rem is in fact a group of closely related results valid for a 
large class of non-equilibrium systems [2l|, [H, [H, [H, HI] ■ 
In a nutshell, FT states that the probability distribution 
for the entropy production rate obeys a symmetry re- 
lation. The Theorem becomes particulary relevant for 
small systems in which fluctuations are large. For this 
reason, the FT has been verified in a number of beauti- 
ful experiments with small systems such as biopolymers 
and colloidal systems ^6|. In the specific case of molec- 
ular motors, FT leads to constraints on the operation of 
molecular motors or nano-machines far from equilibrium 
in a regime where the usual thermodynamic laws do not 
apply [23, Hani. 

In a recent work [s^], we have extended a two-states 
discrete stochastic model introduced in Ref. by in- 
cluding an important variable, namely the number of 
ATP consumed. We have shown that this extended 
model satisfies FT, and we have constructed a thermody- 
namic framework allowing us to characterize quantities 
like the average velocity, the average ATP consumption 
rate and its thermodynamic efficiency. We have also an- 
alyzed the different thermodynamic modes of operation 
of the motor as functions of generalized forces arbitrar- 
ily far from equilibrium. Using FT, we have quantified 
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the "violations" of Einstein and Onsager relations. The 
deviations from Einstein and Onsager relations can be 
studied by considering the linear response theory in the 
vicinity a non-equilibrium steady state rather than near 
an equilibrium steady state. After determining the pa- 
rameters of our model by a fit of single molecule exper- 
iments with kinesin we have formulated a number 
of theoretical predictions for the " violations" of Einstein 
and Onsager relations for this motor. 

In this paper, we further extend the analysis of this 
model. In particular, we provide a more detailed study 
of the modes of operation of the motor and its thermody- 
namic efficiency in relation with the experimental data of 
kinesin. This part contains important information which 
could not be presented in Ref. [30| due to limited space. 
The rest of the paper is devoted to bringing together dif- 
ferent formulations of FT, explaining their connections 
and their physical implications. One way to formulate 
FT involves the generating function of the currents. We 
show that this formulation leads to two versions of FT, a 
long time version of FT similar to the Gallavotti-Cohen 
relation which holds quite generally, and a finite time 
version which is analogous to the Crooks-Evans tran- 
sient Fluctuation Theorem [13, HH] which holds under 
restricted hypotheses for the initial state. Another for- 
mulation of FT takes the form of a property of the large 
deviation function of the current. There are very few non- 
equilibrium models for which the large deviation function 
of the current is known exactly. Our model is sufRciently 
simple for this analytical calculation to be possible, and 
by carrying it out we show that the large deviation func- 
tion of the current satisfies indeed an FT relation. A 
by-product of this calculation is a third formulation of 
FT in terms of the ratio of probabilities for observing a 
velocity v to the same probability for observing a velocity 
—V. The prediction of this ratio of probabilities, which is 
obtained from the large deviation function of the current 
is one of the main results of this paper. We also study 
the connections between FT and the notion of entropy 
production. This entropy production can be explicitly 
calculated using the notion of affinities associated with a 
cyclic evolution of the mechanical and chemical variables, 
and the affinities precisely enter all the formulations of 
FT. We show that the entropy production can also be ob- 
tained from an evaluation of a quantity called the action 
functional [23| . The entropy production and therefore 
also FT depends on the coarse-graining of the description 
which we illustrate by considering three different levels 
of description: purely mechanical, purely chemical and a 
combination of the two. 

The paper is organized as follows: In Section II, we in- 
troduce the model, in Section III we consider the modes 
of operation of the motor, its thermodynamic efficiency 
and the comparison of the model with experimental data 
for kinesin. Section IV is devoted to the formulation of 
FT in terms of generating functions, with its long time 
and its finite time versions. Section V discusses the for- 
mulation in terms of the large deviation function of the 
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FIG. 1: A schematic of the rates for this two-state stochastic 
model of a single processive molecular motor. The position of 
the motor is n and y is the number of ATP consumed. The 
even and odd sites are denoted by a and b, respectively. In 
the case of two headed kinesin, site a represents a state where 
both heads are bound to the filament, whereas site b repre- 
sents a state with only one head bound. Note that the lattice 
of a and b sites extend indefinitely in both directions along 
the n and y axis. All the possible transitions are represented 
with arrows on this particular section of the lattice. 




b b 



FIG. 2: Cycles associated with the evolution of the motor 
of Fig. [T] Left: cycle for the position variable n; the aver- 
age length n undergone by the motor corresponds to half the 
number of full turns of the clock (the factor half is due to the 
period being equal to two lattice units) and with the rates as 
shown. Right: cycle for the chemical variable y, the average 
number of ATP units corresponds to the number of full turns 
with the rates shown. The affinities associated with these cy- 
cles are given by Eq. [TT] for the mechanical variable and by 
Eq. [B5]for the chemical variable. 

current and the significance of the third formulation of 
FT in terms of a ratio of probabilities, and the last sec- 
tion contains the discussion of the entropy production. 



II. A TWO-STATE MODEL FOR MOLECULAR 
MOTORS 

A. Construction of the model 

As a result of conformational changes powered by hy- 
drolysis of ATP, a linear processive motor, like kinesin, 
moves along a one-dimensional substrate (microtubule). 
The state of the molecular motor may be characterized 
by two variables: its position and the number of ATP 
consumed. To model the dynamics, we consider a linear 
discrete lattice, where the motor "hops" from one site 
to neighboring sites, either consuming or producing ATP 
(see Fig.[T|). An alternate representation of the dynamics 
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can be built in terms of cycles (see Fig. [2]). The position 
is denoted by a; = nd, where 2d « 8 nm is the step size 
for kinesin. The even sites (denoted by a) are the low- 
energy state of the motor, whereas the odd sites (denoted 
by b) are its high-energy state; their energy difference is 
AE = ksTe, where fcs is the Boltzmann constant and 
T is the temperature. This model is suitable to describe 
a two-headed kinesin walking on a microtubule, with the 
high-energy state corresponding to the state where a sin- 
gle head is bound to the filament, and the low-energy 
state corresponding to the two heads bound to it. Be- 
cause of the periodicity of the filament, all the even (a) 
sites and all the odd (6) sites are equivalent. 

The dynamics of the motor is governed by a master 
equation for the probability, Pi(n,y,t), that the motor, 
at time t, has consumed y units of ATP and is at site i 
{— a, b) with position n: 

dtPi{n, y,t) ^ - {uJi + Toi) Pi{n,y,t) + 

^ [tJlP,{n + \,y^l,t) + ^lP,{n-l,y^l,t)] ,(1) 

i=-l,0,l 



with i ^ ii 1^ j = 'Ylii ^ j 9-^^ ^ ] = 'Ylii ^ j- We denote 
by tjj and ~ujj the transition rates from a site j to a 
neighboring site to the left or to the right, respectively, 
and with I (= —1,0, 1) ATP molecules consumed. Note 
that the index i is directly linked with the parity of n. We 
expect the rates tJ j and j to be different even if no load 
is applied to the motor, because the interaction between 
the motor and the filament breaks the left-right (spatial) 
symmetry. This requirement is essential for generating 
directed motion. 

Another essential requirement for generating directed 
motion is to break time-reversal symmetry. Such a sym- 
metry is always present at equilibrium, where the de- 
tailed balance condition holds. Detailed balance is bro- 
ken in molecular motors due to chemical transitions in- 
volved in the mechano-transduction process. To obtain a 
simple description of this process, we assume that tran- 
sitions between the two states Ma and M;, of the motor 
[m are possible via two different chemical pathways: 



M„ 



ATP 



M, + ADP 



and 



M 



(2) 



(3) 



where aa (resp. at) are forward (resp. backward) rates. 
The first pathway {a) represents transitions of the motor 
accompanied by ATP hydrolysis, which we call "active", 
and the second pathway (/?) represents transitions driven 
by thermal activation, which we call "passive". In the 
representation of fig. [U the a-pathway represents oblique 
transitions which change both y and n whereas the /3- 
pathway is associated with horizontal transitions which 
change only n. It is straightforward to generalize the 



model with more chemical pathways, but here we focus 
only on these two [ll| . In the absence of an external force, 
transition state theory of chemical reactions requires that 

lallli 



and 



f3a 



(4) 



(5) 



Taking aa = io^^.ab ^ uJi \ = Lot'^^Pb = tob^ and 
"a = £^\afc = ^b~^,l3a = io^°,Pb = we construct 
the transition rates from only four unknown parameters 
a, a' , CO and uj' as follows, 



UJ~b^ 
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= L0', 



(6) 



and with tob^ = ^jj^,^^ — toa^^ — uJb^ — 0. The only ther- 
modynamic force driving the chemical cycle is the free 
energy of hydrolysis. This is quantified by the chemical 
potential AJl = /csTA/i, which is defined by the standard 
expression 34] 



Ap, — fc^Tln 



f[ATP] [AD Pig [P] 



eq 



\ [ATP],q[ADP][P] 



(7) 



where [..] denotes concentration under experimental con- 
ditions and [..]eq denotes equilibrium concentrations. The 
chemical potential of the hydrolysis reaction introduces 
a bias in the dynamics of the motor, which is responsible 
for breaking the time-reversal symmetry associated with 
the detailed balance condition (which holds at equilib- 
rium) . 

Following Ref. [l^, the transition rates can be gen- 
eralized to include an external force Fe according to 
"uJliF,) = t7/(0)e-''r/ and w/(Fe) = w/(0)e+''^/, 
where / = Fed/ (kBT) and 9f are the load distribution 
factors. These load distribution factors take into account 
the fact that the external force may not distribute uni- 
formly among different transitions [s^ . Thus, we may 
write the non-zero rates in the presence of force as: 



LUb 



a e 



UJb" 



LUb 



a e b ■ 



(8) 



In the above expressions, the values of the parameters 
of are arbitrary except for the following constraint: Af- 
ter one period, the work done by F^ on the motor is 
-Fe2d, implying that 0+ + 0^ + 9' + O'^ = 2. Indeed, 
as shown in figlH the simplest model with all the ^^'s 
equal to 1/2, which was studied in Ref. [l3|, does not re- 
produce the experimental curves of velocity versus force 
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for kinesin. The fact that the Of^s are different from 1/2 
agrees with standard models of kinesin, in which several 
chemical transitions are involved, and the force must be 
split unequally among the different transition rates 
We note that this splitting of the force (which involves the 
actual value of the Of^s) is a matter of kinetics, whereas 
thermodynamics enforces only Eqs. HUH The expression 
of the rates given in Eqs.[S][5]is essential for the analysis 
which we develop below: we emphasize that these expres- 
sions are based on first principles. Once these rates are 
decomposed into an active and a passive part, the ratio 
of the passive transition rates in Eq. O follows from the 
condition of micro-reversibility (detailed balance), while 
the ratio of the active transition rates in Eq. [4] requires a 
more general principle for non-equilibrium chemical reac- 
tions. Such a principle is based on the notion of afhnity, 
introduced by de Donder in Ref. 33] to characterize non- 
equilibrium chemical reactions. The de Donder equation 
relates the forward and the backward reactions rates ~a 
and a" of an elementary step, as a consequence of tran- 
sition state theory: 



(9) 



where A is the affinity, defined as — (dGjdQ)T,p in terms 
of the Gibbs free energy G and ^ the extend of reaction. 
At equilibrium, A = and Eq. [5] leads to "a = which 
is the principle of micro-reversibility. Note that Eq. |4] is 
indeed of the form of Eq. [51 with the choice a = aa, 
a = ab and = A/ksT when e = 0. Thus we can 
consider Eq.|3]as a De Donder relation, which generalizes 
the condition of micro-reversibility far from equilibrium 
[33| . Equivalently one can also interpret this equation 
as a particular case of generalized steady state balance 
conditions, we shall come back to this point in section 
IV. 



B. Effective description of the dynamics 

Let us now analyze further the conditions for directed 
motion for this model, which as we mentioned earlier 
are required to break both the spatial symmetry and the 
symmetry associated with detailed balance. These con- 
ditions for directed motion can be derived by construct- 
ing an effective dynamics, which holds at long times and 
large length scales [l3|- Let us first consider a coarse- 
grained description in which the position variable n is 
the only state variable. The chemical variable y may 
not be accessible or we simply do not wish to include it 
in this description. The dynamics of the motor is then 
described formally by a master equation, which can be 
obtained from Eq. [T] by integrating out over all possible 
values of y. We are then left with a coarse-grained master 
equation for Pi{n,t) = J dyPi{n,y,t), which is 

dtPiin,t) = - iuJ., + lJ.,)Pi{n,t) + 

[ tJj P^{n + I, t) + Tu ^ P^{n ^ I, t)] , {10) 



with the same rates as before. Note that these rates 
may still depend on the ATP concentration. As shown 
in Ref. [l3| , an effective potential can be constructed for 
this problem by eliminating one of the sites (a or b) from 
the master equation, Eq. llOi and describing the remain- 
ing dynamics in terms of an effective potential. This 
is the effective potential under which a random walker 
satisfying detailed balance would exhibit the same dy- 
namics. Of course, the same effective evolution equa- 
tion applies to occupation probabilities of site a or b. 
This reasoning pl| gives the effective energy difference 
AE = E{n -I- 2) — E{n) between site n and site n + 2, 
which we write as AE — 2A;bT\E' with 



2 \UJ a^^b 



When the rates of Eq. [5]are used, we find that 



(11) 



,p_ll^J (« + ^) {a'e^^ + ^') 
2 \ {ae^t' + u) {a' + Lo') 



(12) 



Note that the effective potential is independent of the 
load distribution factors and is identical to the ex- 
pression obtained in Ref. [ij except for the change in the 
sign of the force [35]. A nice feature of Eq. [T^lis that the 
conditions for directional motion can now be immediately 
obtained from it, in a way that is completely analogous 
to what is done for the ratchets models in Ref. [l3l- Di- 
rected motion is only possible if the effective potential is 
tilted i.e., AE ^ 0. Thus directed motion requires: (i) 
an asymmetric substrate which means either a ^ a' or 
Lo ^ uj' , and (ii) breaking of the detailed balance condi- 
tion, so that either A/z ^ or / 7^ 0. When A^ = / = 
the system is in equilibrium, the effective potential is flat 
{AE = 0) and no directional motion is possible. A dif- 
ference between this model and with the various ratchet 
models of Ref. ^1^, is that in ratchets the position of 
the motor is a continuous variable. In the classification 
of ratchet models given in (Tlj . our model corresponds 
to a system of class A for which diffusion is not neces- 
sary for motion generation. In this class of models, the 
two ratchet potentials are identical and shifted with re- 
spect to each other in such a way that each chemical 
cycle generates with a high probability a step in the for- 
ward direction. As a consequence of this construction, 
one should expect (and indeed we will find) that in this 
model there is a strong coupling between the chemical 
and mechanical coordinates, and the motor has a strong 
directionality and a large thermodynamic efficiency. 
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III. MODES OF OPERATION OF THE 
MOLECULAR MOTOR, FIT OF 
EXPERIMENTAL CURVES AND 
THERMODYNAMIC EFFICIENCY 

A. Description of the dynamics using generating 
functions 

In this section, we analyze the long time behavior of 
our model using generating functions, which has the ad- 
ditional advantage of making the symmetry of the Fluc- 
tuation Theorem apparent as shown in the next section. 
Let us introduce the generating functions: Fi{X,j,t) = 
X^Ti ^ whose time evolution is gov- 

erned by: dtFi = Mij Fj, where A4[A,7] is the following 
2x2 matrix which can be obtained from the master equa- 
tion of Eq. [TJ 



Vv=0 
"r=0 \ 


- 




\ c ^ 







FIG. 3: Four modes of operation of a molecular motor, as 
delimited by w = and r = [lol |. The lines are generated 
with parameters that we have extracted by fitting the data 
for kinesin in Ref. to our model, and this fit is shown in 

fig-H 



- t^a 



'7b 



Pa 



'LUb - UJb 



(13) 



with ~pn{j) = J2i '''i and p'„(7) = J2i ^'e" 

For t — !■ oo, we find 



{ 



) =^i^,(A,7,t)~exp i§t).. 



(14) 



The method gives also access to higher moments of n{t) 
and y(t). The second moments for instance can be ex- 
pressed in terms of the diffusion matrix 



2A^^ 



[0,0], 



(20) 



where § = i3[A,7] is the largest eigenvalue of M. This 
eigenvalue, d, contains all the steady-state properties of 
the motor and its exact expression is given by: 



with the understanding that zi = A and Z2 ~ 7. These 
first and second moments can also be obtained by cal- 
culating the average of n{t) and y{t) directly from the 
master equation Eq. [T] 1^, |lj, i36|] . 



(15) 



with the notations uja = io a + ^ a and tob = ^ b + ^ b- 

The average (normalized) velocity v is the current of 
the mechanical variable, which is given by 



V — lim 

t — >oo 



< n{t) > 



t 



(16) 



and similarly the average ATP consumption rate r is the 
current of the chemical variable, which is given by 



lim 

t — *oo 



< vit) > 



(17) 



B. Modes of operation of the motor 

From the conditions of vanishing of the currents: v — 
and r = 0, we can construct a full operation diagram of 
a motor, as shown in Fig. [3] for the case of kinesin. The 
curves v — and r — define implicitly / — /st(A^) 
(the stalling force) and AjU — A/ist(/), respectively. The 
stalling force is 



1 ( ja + u;) {a'c^f^+cj') ' 
/st(A^ = - In - — ^ / 

2 \ (ae"^ + uj) (a' + lo') 



(21) 



From Eq. [TH we see that v = —d\d[Q,Q] and r 
— 9-yZ?[0,0], and from Eq. [T5]we find explicitly that 



b — 1^ b 
TOa + 'tOb + '^a+^b' 



UJb 



(18) 



(19) 



which means that for / = /st(A^), v = Q and = 0, 
where ^ was defined in Eq.[T21 At the stalling point, the 
mechanical variable is equilibrated but not the chemical 
variable. Therefore, in general, the motor consumes ATP 
i.e., r 7^ 0, even if it is stalled (in fact, it is only at 
equilibrium, / = A/i = 0, that both v and r vanish). 
Near stalling for A/i 7^ 0, the motor evolves in a quasi- 
static manner but irreversibly. A similar phenomenon 
occurs in thermal ratchets [33, [111 . 

Likewise, the condition A/i = A/ist(/) means that r = 



6 



0. The exphcit form of A^st(/) is 

A^istif) = In -r^ V. 

(22) 

The four different regimes of operation of the motor, dis- 
cussed in Refs. IS] for ratchet models can be recov- 
ered here. In Region A, where rA/i > and fv < 0, the 
motor uses chemical energy of ATP to perform mechani- 
cal work. This can be understood by considering a point 
on the y-axis of Fig. [3] with A/i > 0. There we expect 
that the motor drifts to the right with v > 0. Now in 
the presence of a small load / < on the motor, we ex- 
pect that the motor is still going in the same direction 
although the drift is uphill and thus work is performed 
by the motor at a rate W = —fv > 0. This holds as 
long as / is smaller than the stalling force, which defines 
the other boundary of region A. Similarly, in Region B, 
where rAfj, < and fv > 0, the motor produces ATP 
from mechanical work. In Region C, where rA/i > and 
fv < 0, the motor uses ADP to perform mechanical work. 
In Region D, where rA/i < and fv > 0, the motor pro- 
duces ADP from mechanical work. It is interesting to 
note that the large asymmetry between regions A and C 
in Fig. [3] reflects the fact that kinesin is a unidirectional 
motor. Furthermore the regions A and B do not touch 
except at the origin. With kinesin operating in normal 
conditions in region A with A/i ~ 15, the presence of a 
gap between regions A and B means that kinesin should 
not be able to switch into an ATP producing unit (region 
B), and indeed this has never been observed experimen- 
tally. Note that the explicit expressions for fst and Afigt 
obtained in this model do not depend on the energy dif- 
ference e between the two states, due to a cancellation of 
the numerator and denominator in Eqs. 11811191 Thus the 
diagram of operation of the motor is valid for arbitrary 
value of e. 



C. Fit of experimental curves of velocity versus 
force for kinesin 

We now discuss how the parameters of the model were 
determined using experimental data obtained for kinesin. 
In fig. m (which is also fig. 4 of Ref. [l^l), we have fitted 
velocity vs. force curves for two values of ATP concen- 
trations, and also several curves of velocity vs. ATP con- 
centration at different forces using the data of Ref. 
We have assumed that e"^^ = ko [ATP], which is well 
verified at moderate or high levels of ATP. At low con- 
centration of ATP, there is no such simple correspon- 
dence because it is no more legitimate to treat the ADP 
and P concentrations as constant. We think that this 
is probably the reason why the fit is not as good for 
the lowest ATP values (this concerns the first point in 
the curve at = —5.63 pN and low ATP value in 
fig. m. Nevertheless, we can fit very well the majority 



of the experimental points with this model and we ob- 
tained the following values for the parameters: e — 10.81, 
ko = 1.4 • 10VM"\ a = 0.57 s-i, a' = 1.3 • 10"^-!, 
UJ = 3.5 s-i, uj' = 108.15 s-i, 0+ = 0.25, 0- = 1.83, 
0^ — 0.08, and 0^ — —0.16. These values are reasonable 
within the present accepted picture of the nano-operation 
of kinesin ^8*1. Indeed, e and fcj^^ represent, respectively, 
the typical binding energy {'^ 10 ksT) of kinesin with 
microtubules and the ATP concentration at equilibrium 
(~ 10"^/LtM). Moreover, 0- = 1.83 indicates that the 
back-steps (transitions a — > 6) of kinesin contain most 
of the force sensitivity @]. Furthermore, our framework 
allows us to estimate a maximum stalling force of — 7pN. 

A useful quantity to consider is the distance £, which 
the motor walk using the hydrolysis of one ATP molecule, 
which is 

t^-. (23) 
r 

We find for this model in agreement with Ref. that 
t ~ 0.97(2(i) in the absence of load, which corresponds 
to one step (8nm) per hydrolysis of one ATP molecule. 
Thus the coupling ratio of kinesin is indeed independent 
of ATP concentration and is 1:1 at negligible loads. We 
also find a global ATP consumption rate of r ~ 111 s^^, 
in excellent agreement with known values Q. It should 
be remarked that the global ATP consumption rate mea- 
surements done in ATPase assays in the bulk are in agree- 
ment with the single molecule experiments, which are 
intrinsically very different experiments. Kinesin is well 
described by tightly coupled models which incorporate a 
single mechanically sensitive rate and this is consistent 
with our findings that there is only one transition (tran- 
sitions a — + 6) that has all the force sensitivity i.e., the 
largest load distribution factor. In principle, by changing 
the parameters of the model, we could characterize mo- 
tors which are less tightly coupled, but there will always 
be some coupling because, by construction, the mechan- 
ical steps are intrinsically linked with the chemical cycle. 

We have compared our fit with that carried out by 
Fisher et al. in Ref. [l^] where the same data was fitted, 
and we observe that the outcome of both fits is compara- 
ble. In this comparison, there is an issue of complexity of 
the model under consideration to be taken into account. 
This is especially important in fitting experimental data 
of kinetics, which is typically hard to fit because one has 
many parameters to fit in an expression which is a sum of 
exponential functions. The model of Ref. [TJ] is of higher 
complexity because it involves 4 states instead of 2 states 
for our model, thus we might be tempted to say that our 
model does better in fitting the same data with less com- 
plexity. We believe that this is true when considering the 
data for the velocity only, but if we were to include also 
the data for the diffusion coefficient (which is related to 
the randomness parameter defined in Ref. p^). we agree 
with Ref. that a model with 4 states would then do 
better than a model with only 2 states. 
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FIG. 4: Kinesin velocity vs. ATP concentration under an ex- 
ternal force [13 . The solid curves are the fits of our model 
to data from Ref. From the top down, the plots are 

for = —1.05,-3.59, and — 5.63pN, respectively. Inset: 
Kinesin velocity vs. force under a fixed ATP concentration. 
The solid curves are fits to the data of Ref. From 
the top down, the plots are for [ATP] = 2mM and 5 /iM. 
From this fit, we obtained the following parameters for our 
model : e = 10.81, ko = 1.4 ■ ICp ^W'^ , a = 0.57 s-\ 
a' = 1.3 • lO^'^s^V tj = 3.5s~\ tj' = 108.15s"\ 6+ = 0.25, 



1.83, 



: 0.08, and 9: 



-0.16. 



D. Thermodynamic Efficiency 

Another important quantity that characterizes the 
working of a motor is its efficiency [Tol . Wf\ . In region 
A, it is defined as the ratio of the worls performed to the 
chemical energy: 



Afi' 



(24) 



By definition, rj vanishes at / = and at the stalling force 
/st(A/i). Therefore, it has a local maximum ?7m(A/i) for 
some fraiAfj) between fst < fm < 0. Near equilibrium, 
ryni(A/z) has a constant value, along a straight line 
/m(A/i) c>c A/i inside region A [lC|. However, far from 
equilibrium, the picture is drastically different. We find 
that (i) /rn(A^) is no longer a straight line, (ii) ?7m — ?7m °^ 
A/z for small A^, and (iii) ij^ ~ l/A/.t for large A/i. 
Therefore, rym must have an absolute maximum at some 
Afi > 1. One can also consider the curves of equal value 
of the efficiency within region A. In the particular case 
of the linear regime close to equilibrium, these curves are 
straight lines going through the origin [lO| , but in general 
far from equilibrium these curves are not straight lines 
as can be seen in Fig. [5l and the maximum efficiency is 
reached at a point within region A. 

Note that r/m is substantially larger than r]^. For in- 
stance with the parameters used in fig. [51 the maximum 
efficiency is around 0.59 while rj^ ~ 0.03 (see also Fig. 3b 
of Ref. [201 which contains a plot of rym as a function of A/i 
under the same conditions). Hence, this motor achieves 
a higher efficiency in the far-from-equilibrium regime as 
was also found in other studies of molecular motors us- 
ing continuous ratchet models (see e.g., Ref. [HI). Under 



FIG. 5: Curves of equal efficiency rj within region A (which is 
delimited by the solid line and by the y axis) . The parameters 
are those used in fig.|3]and obtained from the fit of fig. [31 From 
the outside to the inside the curves correspond to r; = 0.2, 
T] = 0.3, r? = 0.4, r] = 0.5 and rj = 0.58. The absolute 
maximum efficiency for these parameters is about 59% and is 
located at A// ~ 14 and / ~ —4.9. 



typical physiological conditions (A/i ~ 10 — 25 fc^T), ki- 
nesin operates at an efficiency in the range of 40 — 60%, 
in agreement with experiments It is interesting to 
note that kinesins operate most efficiently in an energy 
scale corresponding to the energy available from ATP 
hydrolysis. 



IV. FINITE TIME AND LONG TIME 
FLUCTUATION THEOREM 

A. Long time FT 

We note that the rates of Eq. [5| satisfy the following 
generalized detailed balance conditions: 



(25) 



^/poqg-(e:+e-)/-AM/^ (26) 



for 1 = 0,1. Here, P^^^ = l/(l + e-"=) and P^^^ = e-'/{l + 
e^*^) are the equilibrium probabilities corresponding to 
/ = and Afi = 0. We note that these relations, Eq. [^D 
and Eq. 1261 while valid arbitrarily far from equilibrium, 
still refer to the equilibrium state via the probabilities 
P/^'. Using the definition of the equilibrium probabilities, 
one can in fact rewrite Eqs. [^511^ as 



In 



U)b~ 



In 



UJb 



-I' 



OJa 



9+)f~Atil, (27) 
e-)f-Afil', (28) 



for 1,1' = 0,1. Note that these relations are analogous 
to the De Bonder relation of Eq. [5| and to the transi- 
tion state theory equations of Eqs. [IHSl Moreover, by 
combining these two equations, using the constraint that 
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the sum of the load distribution factors is two and then 
multiplying the result by kgT, one obtains 

kBTln'^^j^=F,{2d)-A^{l-l'), (29) 

UJaWb 

which has the form of the steady state balance condition 
discussed in Refs [? ? ]. As pointed out in these ref- 
erences, by identifying the left hand side of Eq. with 
the heat delivered to the medium, i.e. with the change of 
entropy of the medium, the right hand side of Eq. 1291 can 
be interpreted as the sum of the mechanical work Fe{2d) 
and the chemical work — A/i (l — I') on that particular 
set of cyclic transitions {1,1')- In that sense, Egs. [7711^ 
can be understood as formulations of the first law at the 
level of elementary transitions. It is interesting to see 
that these steady state balance relations also lead to a 
FT as we now show below. 

Using Eqs. l|^ and P^ . it can be shown that M 
and A4\ the adjoint of A4, are related by a similarity 
transformation: 

MHf-X,Afi--f] = QM[X,j]Q~\ (30) 

where Q is the following diagonal matrix: 



Q = 









peqg(e-+0+)//2 



(31) 



This similarity relation implies that [A, 7] and t [ / — 
A, A/i— 7] have the same spectra of eigenvalues and there- 
fore 



(32) 



which is one form of FT. Since this relation holds at long 
times irrespective of the initial state, it is a Gallavotti- 
Cohen relation [2l|. Such a symmetry is illustrated 
graphically on Fig.[5]for a simplified case where the chem- 
ical variable is absent. 



B. Implications of FT in the linear regime 

Here, we discuss the implications of FT in the linear 
regime, which leads to the Einstein and Onsager relations 
near equilibrium. Differentiating Eq. [321 with respect to 
A and 7, we obtain: 



d-d d-d 



97 



(33) 
(34) 



The response and fluctuations of a motor are quantified, 
respectively, by a response matrix Ay and by the diffusion 
matrix Dij defined in Eq. 1201 The physical meanings of 
Ajj are: An = dv/df is the mobility, A22 = dr/dA^ 
is the chemical admittance, and A12 = dv/dAfi and 




FIG. 6: Graphical illustration of the symmetry of the long- 
time Fluctuation Theorem (for simplicity a model without 
chemical variable has been used). The largest eigenvalue -& 
is shown as function of ^ = A// for different values of the 
normalized force /. The symmetry of the long time Fluctua- 
tion Theorem corresponds to the symmetry of this curve with 
respect to ^ = 1/2. 



A21 = dr/df are the Onsager coefficients that quantify 
the mechanochemical couplings of the motor. Near equi- 
librium, where / and Afi are small, a Taylor expansion 
of the r.h.s. of Eg. [551 and Eq. [SHwith respect to / and 
Afj, leads to 



dX 

d-y 



dX 



0.0 9-fdX 



A/i 



Afj. 



dXd-1 



927 



(35) 



(36) 



Using the definitions of v and r from Eqs. [551[M1 one 
obtains directly 



A?i / + A?2 Am, 
Xl, f + Xl^ A/i, 



(37) 



-Dij, which arc the Einstein 



with X% = d,^d,^§[Q,Q\/2 
relations, and A?2 = d-^dx-d[{),{)\l2 = dxd^d[Q,(S\/2 = 
A21, which is the Onsager relation. Thus, FT describes 
the response and fluctuations near equilibrium ^2^, 28]. 

It is interesting to investigate how Einstein or On- 
sager relations are broken in non-equilibrium situations. 
The " violations" of Einstein and Onsager relations when 
linear response theory is used in the vicinity of a non- 
equilibrium state rather than near an equilibrium state 
were studied in Ref. [s^l- There, we quantified the vi- 
olations of Einstein and Onsager relations, respectively, 
by four temperature- like parameters, , and by the dif- 
ference of the mechanochemical coupling coefficients, AA. 
Of course, these effective temperatures are not thermody- 
namic temperatures: they are merely one of the ways to 
quantify deviations from Einstein relations; similarly our 
definition of AA is just one of the possible ways to study 
the "violations" of Onsager relations: strictly speaking 
there are no real violations since Einstein and Onsager 
relations apply only to systems at equilibrium. We have 
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origin n(0) = y(0) = in an equilibrium state, and we 
calculate the values of the position n{t) and of the chem- 
ical variable y{t) at time t. We denote the initial state 
by the vector 



FIG. 7: Graphical illustration of the symmetry of the finite- 
time Fluctuation Theorem for a simplified model without 
chemical variable. The left hand side of Eq. is shown as 
function of ^ = A// for different times (arbitrary units). The 
symmetry of the finite time Fluctuation Theorem amounts to 
the symmetry of this curve with respect to ^ = 1/2. 



shown in f30'| some of the possible behaviors of Ty and 
AA for a kinesin motor using the parameters of the fit 
discussed above: in particular, we found that for kinesin 
the maximum value of AA is AA ^ 45 pN~^s~^, and that 
at large A^, AA -10pN"^s~^ We also found that (i) 
one of the Einstein relations holds near stalling (a point 
which we justify more precisely in the next section in 
Eq. I59p . (ii) the degree by which the Onsager symmetry 
is broken (AA 7^ 0) is largely determined by the underly- 
ing asymmetry of the substrate, (iii) only two "effective" 
temperatures characterize the fluctuations of tightly cou- 
pled motors, (iv) kinesin's maximum efficiency and the 
maximum violation of Onsager symmetry occur roughly 
at the same energy scale, corresponding to that of ATP 
hydrolysis (~ 20 ksT) jO]. Experimental and theoreti- 
cal violations of the Fluctuation-Dissipation relation have 
been observed and studied in many active biological sys- 
tems [H, H, 0, S 1 but to our knowledge no experiments 
testing the Fluctuation-Dissipation or the Onsager rela- 
tions have been carried out at the single motor level. 



C. Finite time FT 

The similarity transformation pop implies that all the 
eigenvalues of M [A, 7] and A^^ [ / — A, A^ — 7] are identi- 
cal, not just the largest one. This more general property 
allows us to prove a transient FT, because the dynamics 
of the model at finite time involves all the eigenvalues 
of A4 and not just the largest one. The price to pay to 
have a FT relation valid at finite time is that the initial 
state can no more be arbitrary. We show here that the 
relation still holds, in the particular case when the ini- 
tial state is prepared to be in an equilibrium state (which 
corresponds to the condition / = A/i = 0) and when, 
in addition, a specific condition on the load distribution 
factors is obeyed. To see how this comes about in this 
model, we assume that the motor at time t = is at the 



|Fo) = \Fit = 0)> 



peg 

peg 



With (0| = (1,1), the initial state vector is normalized 
since (0|i^o) = 1- Let us introduce Zi(A, 7, t) = e^I^'Tl*, 
the evolution operator for the generating functions Fi. 
By taking the exponential of Eq. [301 one finds that this 
operator also obeys an FT relation 

U^f-X,Afi~j,t] = QU[X,j,t]Q-\ (38) 

We calculate the following average, similar to Eq. [141 

(e"(/-^)»(*)-(AM-7)y(t)) = (0|Z^(/- A,A^-7,t)|Fo), 

= {Fo\U\f-X,Aii-j,t)\0), 

= (Fo|QW(A,7,t)Q-'|0), 

= {0\U{f,^,t)\Fo), (39) 

We have used Eq.[38]to derive the third equality, and the 
final equation requires the condition: {Fo\Q = (0|, which 
is equivalent to Q^^|0) = |^o) since Q is diagonal. Using 
Eq. [22 we find that this relation holds if the initial state 
is in equilibrium and if the following condition holds 

0a+0t = 0++9^. (40) 

Equation [39] is analogous to the Evans transient time 
Fluctuation Theorem [2J| and to the Crooks relation [2^ . 
An important point here is that the initial state must be 
an equilibrium state while the final state at time t does 
not have to be (and in general is not) an equilibrium 
state. Crooks relation can be derived using a path rep- 
resentation of the ratio of forward to backward probabil- 
ities according to a specific protocol, assuming a Markov 
process and using a generalized detailed balance relation 
between successive states. In our case, the equivalent of 
the generalized " local" detailed balance condition needed 
for the proof is Eq. [38] The symmetry of the transient 
Fluctuation Theorem is illustrated graphically in Fig. [7] 
for a simplified case where the chemical variable is absent 
(or integrated out). 



FLUCTUATION THEOREM FOR THE 
LARGE DEVIATION FUNCTION 



Explicit calculation of the large deviation 
function of the current 



Here, we again take advantage of our knowledge of the 
function i9, which contains all the information about the 
long time dynamical properties of the model, to obtain 
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an explicit expression for the large deviation function of 
the current. To simplify the presentation, we consider 
the simplified description, given in Eq. [TOl in which the 
chemical variable y is not taken into account. In this 
case, the generating function is defined by Fi{X,t) = 
J2y J2n e~^"^Pi{n, y, t), and the matrix M[\] becomes 



M[\] 



-uJb - 



(41) 



By definition, d is the largest eigenvalue of A^(A), so 
similarly to Eq. [15] we have 

d{X)= i {-UJa-UJb+iiLOa-^b)'' (42) 

with the notations uJa — ~uj a + "oJ a and tjf, — ~ujh + *&6. 

We have already seen that i9(A) has the property that 
< e"^" >~ e''^^)* for large t. On the other hand, the 
large deviation function G{v) is defined for large time t 
by 



-G{v)t 



(43) 



in terms of 'P{n/t ~ v) the probability to observe a cur- 
rent V after the motor has gone a distance n from the 
origin in a time t. The relation between t?(A) and G{v) 
is 



Xvt 



(44) 
(45) 
(46) 



Using the saddle point method, we find that i?(A) — 
max„(— G(d) — Xv) and thus d{X) and G{v) are Legen- 
dre transform of each other. We have also —G{v) = 
max(z?(A) + Aw), which can be written in parametric form 



= J tdvVij = v)e~ 



d-d 

^(A*) + A*f = ~G{v) 



(47) 
(48) 



Using Eqs. [13^ we find the following expressions for 
G{v) (see Appendix A for details of the derivation) : for 

i; > 



Giv) 



and for w < 0, 



2v \ ^ ' Y-{v) 



Giv) 



Y+{v) 



-vX (v), 
(49) 

-vX+{v), 
(50) 



200 
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100 
V 
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FIG. 8: Large deviation function G{v), the solid line is the 
exact expression using Eqs. I49I50I and the points are the nu- 
merical evaluation of the Legendre transform using Eq. 1481 
For the values of the rates used here, the average velocity, as 
given by Eq. 1181 is u ~ 80: thus, the system is far from equi- 
librium. Note that G{v) is minimum at v and that G{v) — 0. 



where 

r±(v) 

X^iv) 



\(Z{V)±./ZW 



-4 



2" 



In 



Z{v) ± ^W? 



and 



1/2 



,(51) 



(52) 



with the following parameters 

n = AlS al3 b'uj oXj b, (53) 

= [uJa+l^bf -^{^a^b + ^a~^b) , (54) 

and 'J is the effective potential defined in Eq. [TT] As 
shown in fig.[51 the function G{v) has a single minimum at 
the average velocity v = which was defined in Eq. I18[ 
and at this point G{v = w) = 0, which can be deduced 
from Eq.USl Remarkably, although G{v) is a complicated 
non-linear function of v, the difference G{v) — G{—v) is 
a simple linear function of v as required by the Fluctua- 
tion Theorem. Using the expressions and it is 
straightforward to verify that 



G{v) - Gi-v) = 



(55) 



which in turns implies that the ratio of the probabilities 
to observe v or ~v for large t must obey: 



(56) 



From Eq.[55l and the fact that G{v) and i?(A) are related 
by a Legendre transform, we obtain a third formulation 
of the Fluctuation Theorem: 



d{X) =d{~-^ ~ X). 



(57) 
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Near equilibrium, the large deviation function is well 
approximated by a half parabola when both v when v 
are either positive or negative. For -D > 0, this part of 
the large deviation function becomes flatter and flatter, 
when going away from equilibrium (i.e., for increasing 
entropy production). As a result, the remaining part of 
the large deviation function for u < must be linear 
G(v < 0) ~ — ^'w, so that Eq. [55]is obeyed. This hnear 
part for u < and the half parabola for w > can be 
seen in fig. [S] 

It is interesting to note the central role played by 
the quantity 4" defined initially as an effective poten- 
tial, and which now enters the three formulations of the 
Fluctuation Theorem in Eqs. [5M571 Note that these 
equations were obtained for arbitrary forms of the rates 
~uJa,~i^b,'^a and aJf). If we make a specific choice for 
these rates as in Eq. [51 with no chemistry i.e., for 
A/i = 0, we recover ^ = — /, and then Eq . [571 reduces to 
■(?(/ — A) = '^(A), which is indeed compatible with Eq.l5^ 
when there is no chemical variable and no dependance 
on the rates on chemistry. If the rates are those of Eq. [51 
for A/i ^ 0, we obtain using Eq. [T^ the expression 



J.f-fst(At.))vt^ (fort->oo) (58) 



with the stalling force defined in Eq. [511 and related to 
^-by^- = + (A/i). 

Note that an Einstein relation can be obtained near 
stalling, by performing a Taylor expansion of the r.h.s 
of Eq. [571 with respect to — ^, in way similar to what 
was done in Eqs. [5511551 for the derivation of the Einstein 
and Onsager relations. This procedure means that for 



(59) 



which shows that near the stalling force, the Einstein re- 
lation holds in this description where the chemical vari- 
able is absent [301 . 



B. Discussion 

Note that Eq. [58l puts a constraint on the ratio of the 
probabilities of observing a velocity v to the probability 
of observing a velocity —v. These velocities should be 
estimated from the ratio n/t based on an observation of 
the motor running a distance n (or a distance —n), after 
a time t. This relation has been proven here in the limit 
of long time i, but we expect that such a relation will also 
hold at finite time t under some conditions, as suggested 
by our derivation of the transient FT of Eq. [39l Such a 
relation at finite time was also investigated in Ref. [2^. 

Single molecule experiments on kinesin in which back- 
ward steps were studied were performed in Refs. [l^, [3l| . 
In particular it was shown in these references that ATP 
binding was necessary for backward steps, and that the 



ratio of the overall probability of making one forward step 
(whatever the time) to the overall probability of making 
one backward step (whatever the time) is an exponen- 
tial function of the load, which approaches one near the 
stalling force. It is important to point out that this ra- 
tio which was measured experimentally is not the same 
quantity as the left hand side of Eq. [58l although both 
quantities should be related. In view of this, Eq. [551 
should be considered as a prediction for the behavior 
of single motors like kinesin, which to our knowledge 
has not yet been tested experimentally. This suggests 
that it would be very interesting to probe Eq. [551 experi- 
mentally, by trying to compute directly the distributions 
V{j = ±v) for various times. At the same time, it would 
be also useful to study more extensively the behavior of 
motors of various types near stalling as function of the 
ATP concentration. No notable difference could be mea- 
sured in the stalling force at an ATP concentration ImM 
or 10/iM in the experiments of [19] on kinesin, although in 
principle according to general grounds pJ6l ? ] one should 
expect that the behavior of motors near stalling (and in 
particular the stalling force of the motor itself) should 
depend on the ATP concentration and more generally on 
the details of the chemical cycle of ATP hydrolysis. 



C. Other forms of FT relations 

We have seen that the form of the FT relation depends 
on the state variables of the system, or in other words, 
it depends on the level of coarse-graining of the descrip- 
tion. We consider in this paper the following levels of 
description: 

• (I) : The mechanical displacement n is the only state 
variable. 

• (II) : The chemical variable y is the only state vari- 
able. 

• (III): Both variables n and y are taken into account. 

For case (I), the dynamics is described by the simpli- 
fied master equation, Eq. [TOl and Eqs. (|55ll57p are the 
appropriate forms of FT. 

For case (II), we find that FT can be written in the fol- 
lowing forms (we shall omit the large deviation function 
form of FT for cases (II) and (III)): 



which also holds generally for t oo and 



(60) 



(61) 



In the above two equations, — x can be interpreted as 
the affinity 33] associated with a chemical cycle (see the 
representation of the chemical cycle in fig. ^ and next 
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section for a discussion on the notion of affinity). We 
find that x is given by 



hi 



(62) 



where uj„ — 



and tjf, = ~uj{ + tj( for I = —1,0 



or 1. The physical interpretation of x can be clarified by 
using a method similar to [T4| : after integrating out the 
position variable n from Eq. [TJ and decimating over the 
odd or even sites, one can derive an effective evolution 
equation for the occupation probabilities of the remaining 
sites. In this equation, x plays the role of a effective 
potential for the chemical variable. When the rates of 
Eq. [8]are used, we find that this quantity is given by 



X = -A^ + A/i,t(/), 



(63) 



As expected, the conditions for which x vanishes are the 
same as those for which the chemical current r, given in 
Eq. [ini vanishes. 

For case (III), the FT can be written as follows 



-(*t)+Xr)t 



for t — !■ oo and 



-r) 



i?(A,7) = ^?(-*-A,-x-7)- 



(64) 



(65) 



Here, the affinities associated with the mechanical and 
chemical variables are given, respectively, by — ^ and — x. 
Note that these quantities are in general not the same 
as the ones calculated above in cases (I) and (II) (i.e., 
^ 7^ 5* and X 7^ x)- When the rates of Eq. [8] are used, 
^ = —f and X — — A/i, so that Eq. [32] is recovered from 
Eq.[65l 



VI. FLUCTUATION THEOREM AND 
ENTROPY PRODUCTION 

In this section, we discuss the connections between the 
Fluctuation Theorem described in the last section and 
the entropy production [2^. In particular, we show by 
an explicit calculation, that the parameters \E', Xi / and 
A/i that appear in the symmetry relations Eqs. I55II58I 
are identical to the affinities associated with the various 
macroscopic currents (mechanical and chemical) flowing 
in the system Affinities, introduced a long time 

ago in chemical thermodynamics [1^ , represent intrinsic 
quantities that depend only on the microscopic transi- 
tion rates of the system. Thus, the fact that these quan- 
tities also appear in the Fluctuations Theorems, valid 
far from equilibrium, shows a remarkable connection be- 
tween classical thermodynamics and non-equilibrium sta- 
tistical mechanics. 

We shall first discuss the simplified model, in which 
the chemical variable y is not taken into account in the 



description as a state variable (case (I)). The mechanical 
entropy Syiii) then only contains contribution from the 
disorder in the distribution of the mechanical variable n 
and is defined as 



J2 ^P^(n,i)lnP»(n,t) 



(66) 



in units where ks = 1- Using the master equation, 
Eq. [ini one can calculate the variation of SM{t) with 
time: 



dSu 
dt 



(67) 



^^(^,P,(n,t) 

i^j n 



tjjPjin + l,i))ln 



Pj{n + l,t) ' 



where i,j take the two possible values a and b but are 
different from each other. By transforming the last term 
in this equation as follows 



In 



Pjin+l,t) 



In 



~UJiPi{7l,t)) 

^■Pj(n+l,t) 



In- 



the time derivative of the entropy can be rewritten as 
the difference of an entropy production (which is always 
positive) and an entropy flux. Since we are interested in 
a stationary state where dS/dt = 0, both contributions 
must be equal. From such a calculation, one finds that 
the entropy production and the entropy flux in the long 
time limit are given by 



UJ al-^ b — 1-0 b 



In 



LUb 



UJb 



(68) 

where ^ is the effective potential defined in Eq. [TTl and v 
is defined in Eq. [181 According to the general definition 
[28t , we deduce from Eq. [521 that the mechanical affinity 
of the displacement variable is — 

It is interesting to recall that this result can also be 
derived in a different way : in [2^, it was proven that 
the entropy fiux can be calculated by using a fluctuating 
quantity W{t), called the action functional, which can be 
seen as a local measure of the lack of detailed balance on a 
given path at time t. The matrix M{y) that describes the 
evolution of the generating function of (exp(— ^iy(t))) is 
given by 



UJa '^b + ^^b ' 



UJb 



'UJb - UJb 



This matrix is obtained by deforming the original Markov 
matrix by a parameter v. We emphasize that this defor- 
mation is not the same as that used in Eq.UHto calculate 
the large deviation of the currents. The time derivative 
of W(t) is precisely the entropy fiux. Therefore, we have, 
in agreement with Eq. 1681 



n 



M 



dq 
dv 



(0) 



(69) 
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where q{i') is the largest eigenvalue oi J\f{v). Since the 
matrix M has the property A/"^ [v) = J\f{\ — v), its largest 
eigenvalue q{v) satisfies a Fluctuation Theorem 



q{l -v) = q{v) . 



(70) 



Note that Eq. [68] is valid for arbitrary transition rates; 
if we make a specific choice for the rates such as that of 
Eq. El we find that H = (-/ + fst{Afi))v when ^ 0. 
When A/i — 0, we recover 11 — —fv, a well-known result 
[ll| . but our more general expression for A/x ^ shows 
explicitly the dependence of the entropy production on 
a measurable quantity A/i and its connection to the FT 
which we saw in Eq. 1581 

If we now use a description of the model where only 
the chemical variable y is taken into account (case (II)) 
and the total displacement n is integrated out, we can 
define a 'chemical entropy' Sc{t) as follows 



(71) 



i—a,b y 



Calculations similar to the ones described above allow us 
to derive the purely chemical entropy production in the 
stationary state: 



He 



In 



-rx- 



(72) 



the chemical current r and the chemical affinity x were 
defined in Eq. [19] and Eq. [62] respectively. We also note 
that the entropy production in Eq. [75] can be calculated 
using an action functional whose generating function is 
the largest eigenvalue of the Markov matrix suitably de- 
formed m. 

Finally, we can use the complete description of Eq.]!] in 
which both the displacement n and the chemical variable 
y are taken into account (case (III)). In this case, the 
entropy is given by 

i—a,b n y 

Again, if we make the specific choice for the rates of Eq.[8] 
we find that the following well-known result ll| is recov- 
ered for the entropy production: 



fv + rAfj, . 



(74) 



This relation makes explicit the fact that / is the affin- 
ity of the mechanical position variable with the current 
V, and that Afi is the affinity of the chemical variable 
with the current r. We note that these affinities are dif- 
ferent from those found above in the purely mechanical 
and in the purely chemical models, which correspond re- 
spectively to Eq. [11] and to Eq. [53] The fact that the 
expression of the entropy (and hence that of the affinity) 
strongly depends on the level of coarse-graining used in 
a given description should not come as a surprise. The 
two affinities / and Afi appear in the Gallavotti-Cohen 



relation Eq. [32] This suggests that one should be able 
to construct an effective potential describing the evolu- 
tion of the motor in a 2 dimensional phase space of n 
and y, and that this potential should be equivalent to 
the potential of mean force discussed in Ref. [39| . 

We have seen here that the FT for the currents and the 
FT for the entropy are closely related; this fact is true for 
a large class of models as explained in Ref. |23i] . However, 
although the FT for the entropy holds generally for any 
markovian dynamics as shown in Ref. [23|, a FT for the 
currents exists only if the dynamics can be decomposed 
into cycles with well defined affinities [28| . This is why 
the periodicity of the motion of the motor along track 
and of the evolution of the chemical variable was a crucial 
assumption in our derivation of the FT for the currents 
but was not used when deriving the FT for the entropy. 



VII. CONCLUSION 

We have studied a discrete stochastic model of a molec- 
ular motor, which is a minimal ratchet model. We made 
contact in this paper between various formulations of FT. 
Through a detailed analysis of a simple model, we have 
brought out some physical implications of FT for molec- 
ular motors in general and for kinesin in particular. One 
important message is that FT puts constraints on the op- 
eration of a molecular motor or nano-machines far from 
equilibrium. Further experimental work and theoreti- 
cal modelling is necessary to check more precisely the 
implications of FT for molecular motors. For instance, 
it would be interesting to study a molecular motor in 
which both the velocity and the average ATP consump- 
tion rate could be measured simultaneously, or if this 
is too difficult study more extensively the behavior of 
motors near the stalling force as function of ATP con- 
centration. This would allow a study of the violations of 
the Fluctuation-Dissipation at the level of a single mo- 
tor, which would lead to much deeper insights into the 
Mechano-transduction mechanism of molecular motors. 

Due to the broad applicability of the ratchet concept in 
biological systems, we believe that the results of this pa- 
per should be of general applicability: the model could 
describe processive molecular motors of various types, 
nano-machines like enzymes performing chemical cycles 
or polymers which are translocated through a pore under 
the action of a force (for instance the force created by an 
electric field applied to a charged polymer). More gen- 
erally, we hope that the present work illustrates the use- 
fulness of statistical physics of non-equilibrium systems 
for the understanding of active systems, and in particular 
biological systems. 
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APPENDIX A: CALCULATION OF THE LARGE 
DEVIATION FUNCTION G(V) 



16^2 



r- — 



Y + — 
Y 



(A5) 



To obtain an exphcit form for G{v), we take the square 

of Eq. Ha 



[\ = X*]] ^ 



(Al) 



Using Eq. 221 the derivative on the left hand side can be 
written as 



dd _ 1 V'{\) 



dX 4y[7(A)' 



(A2) 



with C/(A) = {uJa - Wb)^ + 4(t7be^ + be~ ^) (tJ + 



After performing the change of variable 



Y 



^2A I ^ a^b 



(A3) 



and using the parameters f2 and S introduced in Eqs. 
I54[ we can write 



(A4) 



We deduce that Z = Y + 1/Y satisfies 



IT 



(A6) 



There are two solutions to this equation but since Z > 0, 
only the positive solution must be retained which is 
Eq. [5H To obtain Y in terms of Z{v), one must solve 
another second order equation Y^ — ZY +1 = 0. This 
equation has two positive acceptable solutions, which are 
the two solutions Y'^^v) of Eq. Ell We have Y'^{v) > 
Y-{v) = l/Y+{v). Using Eq.[Sland Eq.[A2l we see that 
Y~^ (v) corresponds to t; < 0. Similarly, Y~{v) corre- 
sponds to V < 0. Once the relation Y = Y{v) is deter- 
mined, it is easily inverted using Eq. IA3I to yield \*{v), 
which is precisely X^{v) in the second equation of Eq. lST] 
The final expression of G{v) is obtained by substituting 
this result into Eq. |48l 
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